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Abstract 

We discuss some aspects of the approach of the mass parameters by means of the simple cranking model. In particular, it 
is well known that the numerical application of this formula is often subject to ambiguities or contradictions. It is found that 
these problems are induced by the presence of two derivatives in the formula. To overcome these problems, we state a useful 
ansatz and we develop a number of simple arguments which tend to justify the removal of these terms. As soon as this is 
done, the formula becomes simpler and easier to interpret. In this respect, it is shown how the shell effects affect the mass 
parameters. A number of numerical tests help us in our conclusions. 
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I. INTRODUCTION 



Under the assumption of the adiabaticity of the nuclear motion (the shape variations are slower than the single- 
particle motion), a collective hamiltonian can be defined as the sum of the kinetic and potential energy of deformation 

\ E E {/^i' •' + ^ ^^1' •' 

i i 

In this paper we consider only the deformation of the nucleus neglecting thus the rotational degrees of freedom. 
The set .,/?„} specifies a set of deformation parameters of the nuclear surface. They constitute the dynamical 
variables of the motion. The functions Dij {/3i, ., /?„} represent the so-called mass parameters or the tensor of inertia, 
and U {Pi, Pn] is the potential energy of deformation. Both of these functions depend on the deformation of the 
nucleus. These two quantities are especially important in the study of the dynamic of the nucleus such as the nuclear 
fission (lifetime estimates) [l| and the collective levels of the nucleus '5] . 

Usually, the quantity ., /?„} can be evaluated in the framework of the constrained Hartree-Fock theory or by 

the phenomenological shell correction method. The mass parameters Di il /3j, ., /?„} are often approximated by the 
cranking formula |3] or in the self consistent approaches by other models 

The present work focuses mainly on the drawbacks of the numerical aspect of the calculation of the mass parameters 
on the basis of the cranking formula. In particular, it is well known that cranking procedure is somewhat delicate and 
in some cases can even lead to singular unphysical values for the mass parameters. One of the goals of the present 
work is precisely to discuss this point and to propose a correction to this formula. First and foremost, to our point of 
view, it seems that the main mistake is the presence in the formula of the derivatives of the gap parameter and the 
Fermi level. Therefore, the main part of this work will be devoted to the looking for the arguments (or proofs) which 
could justify the removing of these terms. On the basis of theoretical as well as numerical approaches, this survey will 
also be very useful to demonstrate the relative importance of the different terms entering into the calculations and 
also to show how the shell effects affect precisely the mass parameters. We realize that we criticize an old and well 
established formula, but it is only " the bad" part of the formula which is called into question, not the formula itself. 
This "bad part" is the main source of the problems of the formula. We are intimately convinced that this criticism 
will be very useful and will bring new insight to this subject. 

II. THE INGLIS-BELYAEV OR CRANKING FORMULA 
A. The Cranking or Inglis formula for the mass parameters 

The mass (or vibrational) parameters are given by the Inglis formula [H , Q : 

Where \0) , \ M) are respectively the ground state and the excited states of the nucleus. The quantities Em,Eo are 
the associated eigenenergies. In the independent-particle model, whenever the state of the nucleus is assumed to be 
a Slater determinant (built on single-particle states of the nucleons), the ground state |0) will be of course the one 
where all the particles occupy the lowest states. The excited states \M) will be approached by the one particle- hole 
configurations. In that case, Eq. ([2]) becomes (see appendix \X\ ) : 

^ ^ (fcl 9 I/) (Z| 9 Ifc) 
A,{/3i,.,/3n}-2?.2 ^ ^^-J _ ' ' ^ (3) 

l>X,k<\ ^' ^'^ 

A denotes the Fermi level. The above expression can be simplified further if the following simple property is used: 
(efc - ei) {k\ d /dp, \l) = {k\ [Hsp, d/dp,] \l) = - {k\ dHsp/dp, \l) for k ^ I 
so that Eq. ^ becomes 

{k\^\l){l\^\k) 
A,{/3i,.,/3„}^2n2 E T (4) 
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where Hsp is the single-particle hamiltonian. The single particle states are given by the Schrodinger equation of the 
independent-particle model. 

Hsp\v)=t,\v) (5) 
B. The Inglis-Belyaev or cranking formula with pairing correlations 

It must be noted that in Eq. ^ the denominator e; — tk vanishes in the case where the Fermi level coincides with 
two or more degenerate levels. This is the major drawback of the formula. It is possible to overcome this difficulty by 
taking into account the pairing correlations. This can be achieved through the BCS approximation by the following 
replacements in Eq. 

i) the ground state \0) by the BCS state \BCS) . 

ii) the excited states \M) by the two quasiparticle excitations states ly) (here we consider only the even-even nuclei). 
This led Belyaev Q to two types of non- vanishing matrix elements (see also the demonstration given in the appendix 
E]): 



_9_ 



{iy,-n\d/dl3,\BCS) ^ {u^v^,+Uf,v^){v\ — \^i) v ^ fJ. (6) 



1 f)n 

W,-,.\d/dl3,\BCS)^ (7) 

api 

The later is due to the changing of the occupations probabilities with the deformation, i.e. with u^, and v^: 

iii) the energy Eq by Ebcs and Em by the energy of the two quasiparticles, i.e., by Ei, + E^ + Ebcs- The 

quasiparticle energy being E^ — \J [e^, — + and the BCS state is defined from the "true" vacuum |0) by: 
\BCS)=nk {uk + Vkatatt)\G). 

The " Belyaev formulation" for the mass parameters (formula (|B3p of the appendix |B]) is then: 

Beside this formula, there is an other more convenient formulation due to Bes Q modified slightly by the authors of 
Ref. [1]. In the latter, the above two types of matrix elements have been cast under a more explicit form (see again 
the demonstrations in appendices ICl and ID|) . 

piTT 0\ OA 

(^, -v\ d/dp, \BCS) = - A (H ^ 1^) + + (e. - A) — (10) 

The final formula of the mass parameters takes now the form (cf. formula (jDl[) of the appendix |D|) : 



where 



The two quantities of the r.h.s of Eq. (ITT|) are in the adopted order, the so-called "non-diagonal" and the "diagonal" 
parts of the mass parameters. The derivatives are contained in the above diagonal term When the derivatives 
dX/dPi, dA/d(3i cancel the diagonal term reduce simply to the diagonal matrix element — dHsp/dPi \ v). In order 
to facilitate the comparison with other papers, we recall that the cranking formula is usually cast under a slightly 
different form: 
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where the quantity enclose all (and only) the derivatives [1] because the product of matrix element of i?J is 
"displaced" to the first term of the r.h.s. of Eq. (jl3D . 

III. SOME PRECISIONS ON THE MICROSCOPIC MODEL 

A. The Schrodinger equation 

In this work the numerical tests are based on a microscopic hamiltonian (represented by H^p in the previous 
formulae). The latter is defined starting from a deformed Woods-Saxon potential. In order to obtain a realistic 
model, a spin-orbit and a Coulomb (protons case) terms are also taken into account in the hamiltonian. 
Since in the hamiltonian the kinetic energy operator does not depend on the deformation and since the deformed 
average potential is the most important term compared to the spin-orbit and coulomb interactions, the derivative 
of the single-particle hamiltonian (appearing in the above cranking formula) will be approximated by the one of the 
average potential as in Ref.[l[: 

V denotes the deformed Wood-Saxon potential. In our work, the deformation of the nuclear surface and hence the 
one of the mean field V is of the quadrupole type and is given by the well known Bohr parameters iP,"f)- These 
quantities are connected to the elongation and the axial asymmetry of the nucleus. In this special case, our notation 
reduces here to these two deformations parameters: 

{/3i,.,/3„} = {/3,7} (15) 

The solution of the eigenvalues problem is obtained as follows: 

(i) The eigenfunctions \v) of this hamiltonian are expanded onto the basis' functions of the three dimensional 
anisotropic harmonic oscillator. The representative matrix of this hamiltonian is then effectively built in this ba- 
sis. 

(ii) This matrix is then diagonalized by using a large basis (^ 16 — 20 major shells of the oscillator basis). 

One obtains in this way, the single-particle energies spectrum {e^} and the set of the components of associated eigen- 
functions on the oscillator bjisis. All the details of our microscopic model and the corresponding FORTRAN 
program have been given in Ref . [9| . 

Knowing the single-particle spectrum, the next step is to find the gap parameter A and the Fermi level A by solving 
the standard BCS equations. Technical details of these calculations can be found from Ref lO]. As soon as A and A 
are known, it becomes easy to deduce the BCS amplitudes v^, associated with the energy level e,^. Thus, at this 
stage, all the quantities such as e^,\v) , A, A, u^, which are necessary in Eq. (1111) are known. 

B. Deformation dependence of A and A 

Due to the fact that the single-particle hamiltonian Hgp depends explicitly on the deformation, its eigenenergies 
and its eigenfunctions will be also explicitly deformation dependent. To keep in mind that the deformation dependence 
comes only from the Schrodinger equation, it is useful to highlight this explicit dependence: 

if,p(/3, 7) y, {f3, 7)) = £.(/?, 7) W; iP, 7)) (16) 

Although being explicit, the above notation is somewhat cumbersome, therefore, for brevity the dependence on the 
deformation of the above quantities is usually omitted. 

Starting from the single-particle spectrum the gap parameters A and the Fermi level A are solved from the BCS 
equations (|T7| and (IT8)) as soon as the single-particle spectrum {e^} is known. 



2 ^ 1 
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(17) 
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(18) 



Of course, the deformation dependence of the eigenenergies e,y(/3) involves the ones of A and A. However, this depen- 
dence is not obtained "straightforwardly" from the single particle model itself, that is, from the Shrodinger equation 
(ITSl) . but from the BCS equations (TTTl) and (jl8p . From this point of view, these two quantities (i.e. A and A) must 
not be linked to the deformation as the eigenenergies and eigenfunctions of Hgp do. 

Nevertheless, usually they are considered "explicitly" deformation dependent through the expression of their deriva- 
tives obtained by the so-called lowest-order expansion [l^, [l[: 



d\ acp + bdp dX acj + bd 



■1 



dA bcfj — adp (9A bc-y — ad^ 



dp a^ + b^ d-f a^ + b^ 



(19) 
(20) 



with 



a^Y.^E-\ & = ^(e,-A)i?;3, (21) 

c, = (H ^ W) E-\ d, = - A) (H ^ W) E-^ (22) 

c, = J2A (H ^ W) E-\ d, = J2{e. - A) (H ^ W) EZ' (23) 

In the following the expression " the derivatives" which will be used many times in the text means simply the both 
derivatives given by Eq. (ITOl) and (PO)) . 

IV. THE PROBLEM OF THE PHASE TRANSITION: THE SINGULARITY IN THE MASS PARAME- 
TERS OR THE PARADOX OF THE CRANKING FORMULA 

The transition between the normal and the superfluid phase affects generally the spherical magic nuclei under 
changing deformation [T^ . As we shall see, this phenomenon causes the most serious problem to the mass parameters. 
For convenience, in the following, we choose to discuss just only one parameter, namely Z?^^(/3,7). This does by no 
means restrict the conclusions of this study. 

It is well known that the BCS equations have non-trivial solutions only above a critical value of the strength G of 
the pairing interaction. The trivial solution corresponds theoretically to the value A = of an unpaired system. 
In this case, the mass parameters given by (ITlT) reduces to the ones of the formula (|4]), i.e. the cranking formula 
of the independent-particle model. Indeed, when A = it is clear that the quasi-particle energies of pT|) become 
Ev = — X\, and the quantities u^, Vy are then either or 1 so that the non-diagonal part of the right hand side of 
this formula reduces to: 



^^.yy{u.v, + u,v.) dH^ dH^ 



i^>A/j-<A 



ie.-e,f' ' dp op 



which is the same formula that the one given by Eq. This implies the important fact that in this limit (A — >■ 0), 
the diagonal part (i.e. the second term) of the right hand side of Eq. ([TT|) must vanish. However in a few cases for 
which A « 0, contrary to all expectations, it happens that this diagonal term leads to very large numerical values near 
some "critical deformation" .This well-known singular behaviour constitutes undoubtedly unphysical and undesirable 
effects. This is the main paradox of the formula. 

Because the diagonal matrix element (z/| dHsp/dPi \ v) are finite and relatively small (cf. to the numerical examples 
given below), it is immediately clear from Eq. (1121) that it is the derivatives of A and A which cause the problem. 
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In this respect, the formulae and (I^Ul) are subject to a major drawback due to the fact that their common 
denominator can vanish. This can be easily explained because on the one hand, a is proportional to A (see Eq. (|2ip ) 
and therefore vanishes with it and because on the other hand b is defined as a "random" sum of postive and negative 
values (see again Eq. (HI])) depending on whether the terms are below or above the Fermi level. In the literature, 
this problem has been reported many times [ll|- 15 1, but no solution has been proposed. The authors of Ref. [l| 



proc 

and [13 claim that for sufficiently large pairing gaps A the total mass parameter is essentially given by the diagonal 
part without the derivatives, whereas those of Ref. [l^ affirm that the role of the derivatives is by no mean small 
in the fission process. Other studies [T^ neglect the derivatives without any justification. Even if in some cases the 
derivatives introduce small differences, in others, this is not true at all. In fact, it is easy to realize that these terms 
involve serious and insoluble problems. Therefore, after intensive numerical calculations, we have led to ask ourselves 
if their presence could be called into question. If so, their removal should be justified. Consequently, the problem 
amounts to find good arguments for that. This constitutes the first step of this work. Without further ado, let us 
examine this point in the next section. 



V. ANSATZ AND CORRECTION OF THE CRANKING FORMULA OF THE MASS PARAMETERS 

To overcome "the paradox" of the previous section, we need to start from an ansatz or assertion which justify the 
removal of the derivatives of A and A. Before that, we want to recall some important points. First, the derivatives of 
A and A come from the derivatives of the probabilities in the formula (|8]). It is then worth to recall their expressions: 

(25) 



(26) 




These quantities depend implicitly on the deformation /3 through the level (/?) of the spectrum and also through A 
and A which themselves depend implicitly on the entire spectrum {eo(/3), ei(/3), ...} via the non linear equations of 
the BCS theory. First, we will proceed from a point of view which somewhat looks like the Virtual works principle. 
The following ansatz contains two points: 

1) Thus, we suppose a "virtual change in deformation" S/3 in such a way that one and only level undergoes change, 
namely the level e^. Therefore the "virtual change" of the probibilities would be in this case. 

This is very logical because the changes of A and A are basically due the change of the deformation, not straight- 
forwardly but through the change of the single-particle levels. In other words we have to write {dX/dek) {dtk/dp) 
instead {dX/d/3) because the changes of A and A are determined from the knowledge of the spectrum {ek} which 
depends itself on the deformation of the nucleus. Thus from our "principle" the "virtual change of the probability", 
i.e. Suk, is due to only one level, namely "the corresponding" level ek- 

2) It is well known from the BCS theory that the quantities Uk, Vk, A, A, are related to each other self-consistently 
so that the change of 6uk would depend on its own change trough A and A. For example in the BCS theory we have 
A = UkVk and 6 A would depend on the changes duk and 6vk which in turn depend on SA. This constitutes an 
illogical consequence. To overcome this problem, we state the following intuitive assertion which avoid such a problem 
and hence the "initial" problem of the derivatives of A and A.: "Between two infinitely close "virtual" deformations 
/3 and /? -|- 5/3, the quantities Uk, Vk, X and A must verify the simple (but not so obvious) following conditions": 



duk dX duk dA 
'dXd^^'dAd^ 
dvk dX dvk dA 
dX dek dA dek 

in such a way that Eq. 1^7^ and (pS)) reduce to: 



(29) 
(30) 
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Suk = 


duk dek 
dek dp 






Svk = 


dvk dek J „ 
dek dp 






duk 


duk dek 






op 


dek dp 






dvk 


_ dvk dek 






op 


dek dp 




duk 


duk dek 


duk dX 


duk dA 


dj3 


dek dfi 


+ 4 
dX d/3 


dA dp 


dvk 


dvk dek 


dvk dX 


dvk dA 


dfi 


dek dp 


H h 

dX dp 


dA dp 



(31) 
(32) 

Implying that: 

riiti. f)ni. f)pt^ 

(33) 
(34) 

up uek up 

Must be employed instead 

rtiii^ f)']ii. rif^i^ fill I. f)\ fill!. f)A 

(35) 
(36) 

In fact a real change of the deformation implies "all the levels of the spectrum". Consequently the "virtual" changes 
Suk and Svk must be applied "successively" to all the levels {uk = ui,U2. ....). We recall the important fact that Eq. 
(l35t is used in appendix |D] to obtain the formula (ITlT) whereas Eq. (l33l) leads to the following more simple result 
(which consists of simply ignoring the derivatives of the "old" formula): 

In this case, it is very important to note that Eq. ([57| reduces to Eq.Q of an unpaired system and removes the 
previous "paradox" of the "old" formula. 

Furthermore, for convenience it is to be noted that it is possible to collect the both terms in the right hand side 
of the above formula, say, the non-diagonal and the diagonal parts in only one term. Indeed, with the identity 
luyVi, = AE^, we note that for ly ^ fi: 

{E, + E^f {2E,f 8£„ 

so that the missing term in the "non-diagonal part" is precisely the second contribution of the right hand side of Eq. 
([57)1 . Consequently Eg. (1571) can be cast under the more compact form: 



A,(/3) =2h'Y.T. (H ^ Ia^> (a^I ^ 1^) (39) 

with 



(40) 



Thus, the only difference between the two versions of the mass parameters Eq. (fTT|) and Eq. (|37l) or Eq. (jJH) is the 
absence of the derivatives dX/dPi and dA/dPi in the second version of the formula. The latter constitutes for us 
the cranking formula without derivatives or the "corrected" formula. In the following c^f^ will be called simply "the 
coefficient", and for v ^ ^, dHgp/dPi will be referred to as "the squared matrix element". 



VI. NUMERICAL EXAMPLES AND OTHER ARGUMENTS IN FAVOUR OF THE "CORRECTED" 
CRANKING FORMULA (I.E., WITHOUT DERIVATIVES) 

A. Comparisons between the two variants of the formula for A ~ 

In the present section we will try to prove that the modified or " corrected" cranking formula for the mass parameters 
is consistent with "reasonable" results. In this respect, we have first looked for a "critical" case where there is no 
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FIG. 1: (Bottom) Neutron contribution to the mass parameters Bpp for the magic nucleus l^^Xegi', The calculations are 
performed within the cranking formula including the derivatives and for the same formula without derivating; Note the quasi 
divergence of the version with derivatives near the deformation j3 = 0.1; The curve without derivatives is characterized by 
three region, the first with very small mass parameters (no BCS solution), the second with strong ( but finite) values with a 
maximum at 13 = 0.2, and from /3 = 0.3 a mean values about 100 ~ 250h^ MeV~^ ; (Top) same as (Bottom) but for "mid-shell's 
neutron" number, i.e., for the nuclei ge^Xero; Note here, the regular behaviour of the pairing gap and hence the one of the 
(Bfiis ~ 50- ISOfi^MeV"^). 



mass parameters 



BCS solution apart from the trivial solution (A — 0) and applied then the both formulae. Such situation is obtained 
generally for a magic neutron or proton number. Here, we use the nuclei 54^X652 for which we have plotted in fig. 
[T]-bottom the neutrons' contribution (N = 82 = magic) to the mass parameters vs the quadrupole deformation /3 in 
the two versions of the formula. 

In these calculations the deformation parameter of the axial asymmetry is fixed to the value 7 = 0° (prolate shape). 
The mention "with" between parentheses signifies that the calculations are performed with the formula PT|) including 
dX/d/3i and dA/d(3i, and by "without" the same ones done with the formula Eq. ([57]) without these terms. From this 
example, it is clear that whenever the gap parameter A becomes small (here A < 0.2 MEV) the two formula give very 
different results for the mass parameters. The major difference is reached for /3 ~ 0.1 (for which A « 0.004) where 
we obtain values about 6.000000 MEVh~^and 25MEVh~^ioT the two versions of the formula of the mass parameter 
Dfjp. This singular behaviour (of course within the version including the derivatives) has obviously nothing to do 
with reasonable physical values for the mass parameters. However, both formula give close, though non-identical, 
results as soon as A > 0.2 MEV. Thus, this simple numerical test confirms that it is the derivatives which appear in 
the formula (|lip that are responsible of the divergence character of the mass parameters for the case where the valid 
BCS solution breaks down. 

More precisely, we have checked that it is the common denominator in Eq. (fTO)) and I^U^ which causes the problem. 
Indeed, in this denominator the quantities a and b defined by (j21|) vanish simultaneously. On this point, we give in 
the following table further details concerning the numerical values of a, b, cp, dp which intervene in the quotients of 
the derivatives in Eq. dUl), and Eq. I^Uli . 

Of course, it is to be noted that numerically A does not vanish rigorously for the trivial solution. This is due to the 
occupations probabilities which are numerically not exactly equal to or 1. As mentioned before, the large values of 
the derivatives are due to the denominator + b^ which is very small (For explanations see section lIVp . As already 
noted, this constitutes the main defect of the formula including the derivatives. On the oher hand, we have also verified 
that there is no crossing levels near the Fermi level for this case (/3 = 0.1). Thus the derivatives can diverges without 
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TABLE I: Some numerical details concerning the calculations of the derivatives given by the formulae p9p and (|20p. 



p 


7 


a{MEV-'^) 


b(MEV-^) 


cp{MEV-^) 


di3{MEV-^) 


A{MEV) 


1${MEV) 


§{MEV) 


0.1 


o 




0.0096 


0.0034 


0.0190 


-25.7134 


0.0044 


2366.3 


-832.6 



any problem of the so-called crossing levels, as it is often claimed [12[. 
B. Single- particle contributions to the "element" Rl, 

For a valid BCS solution (i.e. for A G, or approximately A ^ 0.80 in our examples) the problem of the divergence 
disappears. Indeed, from fig [l]-bottom as soon as /? exceeds 0.3 the gap parameter reaches its normal value and the 
divergence is no more present, the both formulae with and virithout derivatives give close results. 
Unlike the closed shell, the mid-shell regions (neutrons or protons) are characterized by a very regular behaviour of 
the pairing A and hence of the mass parameters. This is the case in fig[l]-top for the neutron number = 70 of the 
nuclei 54^X670. However, as in the preceding example, the both formulae (with and without derivatives) do not give 
the same results). This is entirely due to the derivatives contained in the expression of which is defined from Eq. 
p2p . Indeed, the quantity, is the contribution of three terms which are summarized with their most important 
corresponding coefficient in the table (|ll]). The energy levels are labeled by 1/ according to an increasing order as 

TABLE II: The most important single-particle contributions to the term given by the formula p2p . Note that the derivative 
d\/dp is not so small compared to the diagonal matrix elements of the derivative (of. text). 



P = 0.50, 7 = 0°, A = 0.90 MeV 




coeff. 


first contrib. 


seed contrib. 


third contrib. 


level number v 




{^\ ap- W) 


d\ 

dp. 


(..-A) c'A 
^ dp. 






{MeV) 


(MeV) 


(MeV) 


40 


0.117 


-0.81 


7.66 


1.03 


41 


0.292 


-14.19 


7.66 


0.38 


42 


0.186 


-14.50 


7.66 


-0.75 


43 


0.108 


0.72 


7.66 


-1.07 



already precised and the associated eigenfunctions are noted by In this particular case, we have took = 82 
neutrons. There is thus 41 pairs of particles in time-conjugate states. For this reason the Fermi level is close to the 
41*'' level, i.e., the last filled level. This explains why the most important levels are labeled with indices which are 
close to 41. From this table it is clear that the third contribution is small near the Fermi level, dA/df3i being also 
relatively small for such a value of the pairing gap (A = 0.90 MeV). On the other hand, the contribution dX/d/3i 
which does not depend on the level v is by no means so small (^ 7.66MeV) compared to the matrix element. This 
result contradicts some works in which the derivatives are neglected, claiming that the derivatives are small in the 
case of a valid BCS state. Further checkings have shown that it is dA/d/Si which is responsible of the singularity 
(A « 0) whereas it is d\/d(3i which modifies the results for valid BCS solutions. Thus, in all the cases the derivatives 
modify the results. In fact, in terms of physical arguments, we can say that the derivative dX/d/3i perturbs or "masks" 
the shell effects of the diagonal matrix element (see the next subsection) and this consitutes an additional physical 
argument in favour of the removal of the derivatives. 

C. Matrix elements: comparison between the diagonal and the non-diagonal parts. The shell effects lie in 
the diagonal part 

At this stage, the effects of the "perturbations" involved by the derivatives are known. Therefore, from now, we 
will concentrate only on the formula of the cranking without derivatives. It is interesting to compare the magnitude 
of the contributions of the diagonal and non-diagonal parts given by Eq. p7p . These two contributions are plotted in 
fig. [2] for a magic number of neutrons = 82 (bottom) and for a mid-shell neutron number = 70 (top). In the both 
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FIG. 2: (Bottom and Top) same as fig 1; Here, the figure details the diagonal contribution without derivatives and the non- 
diagonal one; Apart from the spherical case (/? — 0) for N = 70 which constitutes a "special case", the diagonal part is by far 
always more important. 



cases, the diagonal contribution is dominant (about 10 times larger) if we except the region of the spherical shape 
(/3 — 0) for TV = 70 which seems to be very particular. Apart from the noted special spherical case, the explanation of 
these results can be easily understood. Indeed, first, the diagonal matrix elements are always by far more important 
than the other because in the integrand the derivative of the potential are multiplied by a squared wave function, 

H^I->«Hf l->= / \^A^,Pn)f^^^^if^d'r (41) 

whereas for the non-diagonal matrix elements, i.e., for 

H^lM)«Hf / (p:(^,/?,7) ^^^?/-"V M(^:A7)rf^r (42) 

the two radial wave functions are not always in phase, their product </?* (1^, /3, 7)(/7^(7^, /3, 7) has not a definite sign, 
the contributions to the integral will have tendency to cancel each other, or at least to lead to a small values. 
It is clair that the more important coefficients Cj^^ (diagonal and non diagonal) are those which are the closest to the 
Fermi level. Indeed, these coefficients are which given by (see Eq. (HO])), arc maximums when the quasiparticle energy 

— \J (ci, — A)^ -|- reaches its minimal value, i.e. for e^, « « A. The few important non-diagonal coefficients 
are generally multiplied by weak matrix elements. Indeed, for closed levels (1/ « /i, v ^ ji), the radial wave functions 
are often in "opposite" phase (like two " consecutives" wave functions of the one dimensional infinite square well) so 
that the contribution to the integral will be necessarily very small. This explains why the diagonal part is dominant. 
Before closing this paragraph, we would like to add an important remark concerning the shell effects. The latter 
are of course due to quantum mechanics. They involve sudden variations in the mass parameters with the number 
of particles. They appear when a few states contribute essentially to the mass parameters. In fact, for the mass 
parameters, we have a very limited number of diagonal matrix elements which contribute essentially to the mass 
parameters opposed to an admixture of a large number of small non-diagonal matrix elements. Consequently, it is 
clear that the shell effects will be contained in the diagonal part (of course without the derivatives) . 
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D. Magnitude of the different terms 



We have found numerically that when there is not BCS solution the diagonal term (without derivatives) is practically 
equal to zero (~ 10~^ MEV) so that the mass parameters reduce only to the non-diagonal part of the formula. In this 
case the mass parameters are very small. Because the values and the fluctuations of the mass parameters are mostly 
due to the diagonal contribution, our numerical study will concern exclusively this term. In the following, once more, 
we will concentrate obviously on the version of the formula without derivatives. The "simple" diagonal part can be 
cast under the following form: 

Di3i3{diagonal part) = ^ Cy^ {v, \ — \v) (43) 

y "Pi 

with, 

cu. = (hV4) , (44) 

As mentioned before, only few (about two up to five at all) energy levels contribute really to the final value of the mass 
parameters. Indeed, this function is peaked at A and due to the power of the denominator of Eq. (j44l) decreases very 
rapidly so that only few levels give important coefficients. For a valid BCS solution, a relative weak value of A gives 
significant coefficient c^,y.For example, we give in fig. [3] (bottom and top) the coefficients c^i, (black stars) and the 
square of the matrix elements dV/dfii \v) (black squares) for the neutron contributions of the same two isotopes 
as before and the same deformation. First of course, we especially can check that these coefficients are important near 
the Fermi level (41*'* and SS*'' level). We have thus computed the four major neutron contributions of the energy levels 
and we have found 97% and 82% respectively for = 82 and iV = 70 cases. These contributions are not only due 
to the strongest values of the coefficients c^^ but also to the "corresponding" matrix elements (see Eq. (H5)) ).We also 
can see very clearly the shell structure of the squared matrix elements which is the same for these two cases (because 
the deformation and hence the spectrum are the same in the both cases). An other remark is that the maxima of 
the matrix elements are situated at the levels N°15, 25, 41, 65 which correspond to the nucleon numbers 30, 50, 82, 130 
which are magic or nearly magic, but these numbers are obtained for a weak deformed shape of the nucleus, whereas 
in the binding energy they match the spherical shape. Moreover these maxima are immediately preceded or followed 
by small values. Thus the shell structure of the matrix elements appears less regular than the one of the energy levels. 
As we can see, the shell structure is "contained" not only in the level density but also in the matrix elements. Thus 
the mass parameters are characterized by a complex "dual" shell structure. 

At last, as a special case, we must point out the very strong values of the coefficients (i.e., the stars) in Fig [3]-bottom 
involving very large the mass parameters at (3 = 0.2. This very exceptional situation arises in the so called one 
crossing levels (see Ref . ( [TTj ) and ([l2|)) near the Fermi energy. This will be analyzed in the next subsection. 



E. Effects of the crossing levels on the diagonal part 

In Fig. m we have plotted the energy levels for different values of the deformation parameter /? (7 being fixed to 
0). The BCS Fermi level is indicated by a star. This level is situated in the middle of a gap constituted by the last 
filled and the first unfilled levels. As (3 increases, these two levels approach little by little each other. From the first 
(/3 = 0) and up to the fourth deformation (/3 = 0.2), the Fermi level lies always in the middle of a gap corresponding 
to the neutrons magic number N = 82. In these locations, the level density is so small with respect to the pairing 
strength so that it (i.e., the level density) is not able to insure a valid (non-trivial) BCS solution. In these four cases 
we have thus A « 0.0044 MEV and hence no valid BCS solution. Therefore, due to the weakness of A, the diagonal 
term of (|37p does not contribute at all to the mass parameters. In these situations, the mass parameters are very 
small ( ^ 2bh^ MEV^^in our case). However, for the deformation (3 = 0.2, the two levels become very close to the 
Fermi level. It occurs in the immediate vicinity of the Fermi level a crossing levels. This happens inside a gap and 
gives rise to two very strong coefficients. Moreover, the corresponding values of the squared matrix elements are also 
important, giving practically all the contribution to the mass parameters. This phenomenon is characteristic of a 
phase transition (compare with[T]-bottom). The remaining cases {(3 beyond 0.25) can be considered as "normal cases". 
They correspond to values about w 0.1 — 0.2. 
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FIG. 3: Following the same examples as those given in fig (1) and (2), we precise here the single-particle contributions of the 
diagonal- part for the both isotopes; These contributions are made of products of the coefficients (-A') by the square of the 
diagonal matrix elements (■) (see Rel. (|43p '): The arrow indicate the levels which contribute the most to the mass parameters; 
For example in the nucleus 54^X632 (bottom) we have 82 neutrons. Each level supports two neutrons; The las filled level is 
thus the 41*''. Therefore the fermi energy is situated between the 41"* and the 42*'' level; In the vicinity of this energy we see 
the important values of the coefficients which are represented by stars; This is essentially due to a crossing level near the fermi 
energy (see text); The black squares indicate the values of the corresponding squared matrbc elements; Coefficients as well as 
matrix elements follows their proper shell structure; Thus the mass parameters are governed by a "dual" shell structure (see 
text); (Top) analog results for the second isotope 54^X670 with N = 70 neutrons; In the latter case,of course, the fermi level 
will lie between the 35*'' and the 36*'' level; For this case the coefficients are less important. 



F. Some other important precisions 



We have focused our discussion around the mass parameter D/^p^ but similar studies will lead to analog general 
conclusions (seen above) for the two other mass parameters D^^ and Dp^. However, we will make further useful 
remarks. For the quadrupole deformation (/3,7), we write the detailed form for the three mass parameters (naturally, 
once more without the derivatives): 

D,,^2h'J2J2''-^^ Hf Im)' (45) 
D,,^2h'Y.Y.''-^^ Hl^l/^)' (46) 
Dp, = 2h'J2J2 ^-^^ f Im) (mI ^ W) (47) 

with Cj,^ given by Eq. PO)) . 

It is easy to compare these three quantities. First it is to be noted that they have the same coefficient so that the 
only difference will come from the matrix elements. The latter have not a definite sign in such a way that the third 
mass parameters will be the result of the sum of positive as well as negative terms. It will be thus necessary smaller 
than the two first mass parameters for which we have only positive contributions. Moreover, it is found numerically 
that Dpp is generally larger than D„. This can also be easily explained by the fact that the single-particle potential 
is generally more sensitive to the f3 degree of freedom than the 7's one, especially for the region the Barium isotopes. 
In other words for the quadrupole deformations and a Woods-Saxon potential we have, 

dV/dp > dV/d-f (48) 
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FIG. 4; Energy levels for 54^X682 at differents deformation; For small deformations due to the spherical symmetry there is a 
high degeneracy; For this reason it seems that ther are less levels ai P — and approximately for the nearest deformations. 
The fermi level is indicated by a star; Due to the neutron magic number A'' — 82, it is situated in the center of an important 
energy gap (between the 41**" and the 42*'' level); Note that the gap reduces gradually as P increases from the spherical shape; 
For the four first levels there is no BCS solution and the mass parameters are small; For /3 — 0.2, we obtain a crossing level 
which contributes strongly to the mass parameters through the diagonal part; Compare these comments with the results of 
fig([T}bottom). 



Although the present work is devoted exclusively to the mass parameters, we would like to add an important remark 
concerning the moments of inertia. Indeed, these quantities can also be evaluated with the help of the cranking 
approximation including pairing correlations. (j49p : 

Because the latter looks like to the one of the mass parameter (of course without the derivatives), it is tempting to 
conclude that the behaviour of these two quantities will in principle be the same. However, it is well known that unlike 
the mass parameters, the behaviour of the moments of inertia is always found as a smooth function of particle-number. 
It depends mainly on the deformation of the nucleus. This means that contrarily to the mass parameters, there are 
no shell effects in the moments of inertia. From our previous study, it is easy to interpret this difference. 

First, this is due partly to the lower power of the denominator of Eq. (|49p but the main reason is that there is 
no diagonal contribution {1/ = /i) to this sum because the "coefficient" vanishes when v — fi. Remembering that the 
shells effects come almost exclusively from the diagonal part of the total contribution, we can conclude here that they 
will be necessary quasi absent in the moments of inertia. 



VII. COMPARISON BETWEEN THE BOTH FORMULAE THROUGH THE EXPERIMENTAL COLLEC- 
TIVE LEVELS 



We have also solved numerically the generalized Bohr Hamiltonian with the help of the numerical fortran code 
of Ref. [l7| . The six inertial functions have been calculated of course via the two versions of the cranking formula 
and the potential energy of deformation with the Strutinsky method. The two types of calculation are compared to 
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FIG. 5: Collective levels deduced numerically (sec tcxte) from the generalized Bohr Hamiltonian including the seven inertial 
functions (three mass parameters, three moments of inertia, and the potential energy of the deformation). Two typical cases 
with and without phase transition are considered. 



the experimental levels. For this task we chose two Barium isotopes. This choice is justified by the fact that the 
magic nuclei ^^^Bas2 undergoes the phase transition near the spherical shape whereas for the other (non magic) nuclei 
^^^Bo76 there is not the case. 

The formula without derivatives is called as the "new" formula whereas the other is called the "old" formula. The 
experimental low lying collective levels are denoted by "exp". From the figure it is clear that for the phase transition 
(^■^^i?o) the both formulae give very different results. The advantage of the new version of the formula is clear because 
the old version gives a completely absurd result. In the other case the difference is much more less important but the 
theoretical "new" and "old" formulae do not give identical spectra. 



VIII. CONCLUSION 



The cranking formula of the mass parameters is widely used in the study of the dynamics of the nucleus. However, 
its applications are not free from ambiguities and contradictions. The aim of this paper is to convince that these 
problems are inherent to a spurious presence of the so-called derivatives in the formula. Therefore, by means of a 
simple "proved assertion" based on a number of pertinent arguments (see text), it is simply recommended to remove 
these derivatives. It turns out that the cranking formula without derivatives is no more subject to the cited problems 
and becomes simpler to interpret. Let us re-examine some advantages: 

i) The " new" formula reduces in a naturally way to that of the unpaired system when A = whereas the old version 
leads to an inextricable paradox. 

ii) The problem of the unphysical large values of the mass (singularity) and thereby unphysical collective spectra 
disappears. 

iii) The shell effects are physically no more "masked" or "perturbed" by the derivatives. 

iv) The formula becomes more transparent and a number a features have been highlighted in this paper. For 
example, the shell effects are not only connected to the level density as it is usually claimed, but also to the matrix 
elements which have also their proper shell structure. Thus, the shell effects appear to be governed by a "dual" shell 
structure. Therefore, examining the mass parameters only through the level density, as it is usually done, is surely 
insufffcient to explain the results. 
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We hope that this work will help the reader to a better understanding of the mass parameters and will encourage 
further investigations in this direction. 



Appendix A: The cranking formula in the independent-particle model 

The cranking formula for the mass parameters is given by Eq. ([U ) . If we want to apply the cranking formula in 
the framework of the independent-particle model, we assume that the excited states are simply one particle-one hole 
excitations. We proceed to the following replacements: 

\0) — >■ \core) and \M) ^t>\^k<x |core) where \core) = Yln<x |0) is the ground state or the core of the 
independent-particle model. In this expression, |0) is the "real" vacuum and is the creation operator of a real 
particle in the state n. 

In the second quantization formalism, we can write: 
5/aA=E.,^ W\d/dl3.,\f,) a+a^ 

We choose \fj,) as the eigenfunctions of the single particle hamiltonian Hgp = The are then the associated 

eigenenergies. 

The matrix element entering in the cranking formula reads then: 
(M| d /dPi \0) = {core\ a^^^Qi^^ ^^^^ {v\ ^ a+a^ \core) 

It is clear that the only non- vanishing matrix elements imply that: fi = k and v = I. 
(M| d /dPi \0) = {l\ ^ |fc) for A: < A and I > A. with {core\ \core) = 1. 

Assuming that in the independent-particle model the energy of the nucleus amounts simply to the one of its nucleons 
we will have Em — Eq — ei — ek, therefore the cranking formula becomes: 



l>\,k<\ ^' 



Appendix B: The cranking formula with pairing correlations 

The BCS theory gives the ground state and two-quasiparticle excited states for a paired system as follows 
M--\BCS) = Uk>o(.^k + Vka+atk) |0) , 

a+a+\BCS) = a+ a+Y[^^^^^{u}, + Ufea^a^j.) |0) for v and a+at^\BCS) = {-v^ + w^aja^,.) Dfe^i^K + 

VkCtt^-k) |0) where {ui.,Vi.) are the usual probability amplitudes of the pairs of particles in the mutual time-reversed 
states |A;), \—k). 

Starting again from the cranking formula ©, we have to evaluate the matrix elements :{M\ d jdfii \0). Then, we 
will replace the state \0) by the BCS state. The excited states \M) are supposed to be two quasiparticle excitations 
states \M) = a^a^ \BCS). 

Expressing the operator d jdfii in the second quantization formalism and keeping in mind however that the differential 
operator acts not only on the wave functions of the BCS state but also on the occupations probabilities Uk,Vk, we 
have to consider two types of differentiation: 

\ / wave June \ / occup.prob 

We must therefore to evaluate successively two types of matrix elements 



1. Calculation of the first type of matrix elements 



We have in second quantization representation: 

\ ^ / wave June 

Applying this operator on the paired system we find 

\ ^ / wave June 

Two excited quasiparticle states of an even-even nucleus are obtained by breaking one pair of correlated particles so 
that we must take v ^ fi (It is easy to verify that this correspond to the first type of the excited states given above). 
Using the inverse of the Bogoliubov-Valatin transformation: 
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in the previous expression, we find 

(^1 l'^^ {ut^Vi^a+ati^ + .not.two.qp) \BCS) 
We considered here only the states with two quasiparticle (for the even-even nuclei). We calculate then the first type 
of matrix elements: 

h = (M| E.,^^. {y\ 4 Im) u.v,atal^ \BCS) 

The above form of the formula suggests that the excited states must be of the form \M) = a^a^i \BCS) = \k,—l). 
Wc obtain then: 

h = {BCS\ a_ia^ J2,.,^^u Wl '^^ w^^^a'"^"^^ \BCS) 

= T.u,t,^y {v\ 1/^') '-'''^'^A' {BCS\ a_ia,^a+at^ \BCS) 

We use the following usual Fcrniions anticommutation relations: 

Thus the quantity between brakets gives: 

{BCS\ a_ia^a+at^ \BCS) = {BCS\ a_, {5^k - a+a^) at^ \BCS) 
consequently: 

E,..M#;^ ('^l WI l'^^ "'•'''^^ (Slf^Sr^k - 5-^k^-^l) 

= {k\jl^\l)ukVi-{-k\j^^\~l)u.iv^k 

the condition [i ^ v implies k ^ I otherwise the result cancels for k = I. 
Noting that if T is the time-reversal conjugation operator 
(-fc| 41-0 = T+^T 10 = (fc| T-i^T 10 
and assuming that g|- is time-even, we get: 

(-fclo|-l-0 = (fcl4lO 

Moreover, using the well-known property 

U-l = Ul, V-k = -Vk 

we obtain finally: 

/i = (fc,-^|(^) \BCS) = {ukVi+uiVk){k\-^\l){l-5k,i) (Bl) 

V A-'i / wave func 

2. Calculation of the second type of matrix elements 

Differentiating the BCS state, we obtain: 

im) , l^c-s) = E.( tr + wi'^r^-r) n^^K + Vk4at,) |o) 

\ / occup.proo 

This corresponds to the second type of the excited states given above because for the pair in the state r we have a 
linear combination of the vacuum state with a state of two particles in conjugate states. The excited states will be 

necessarilv here, of the following form: 
\M) = aiat^\BCS) = \m,-m) 
We have therefore to calculate: 

h = {BCS\ a_^a^ E.(f|r + m^'r^^-r) Uk^ri^k + Vka+atk) |0) • 
It is simple to show that: 

YlkM^'k + ^^ka+at,) |0) = {ur + Vra+atX^ \BCS) 
then 

h = {BCS\ a_^a^ Y^rim + fe«r «-x)K + Vratat,)-' \BCS) 

The case t leads to the cancellation of ■ The only non vanishing possibility is given by r = m. 

h = {BCS\ a_^a^{^ + f^a+at„)(w„ + v^a+,atj-^ \BCS) 

Using the inverse transformation of Bogoliubov-Valatin by "selecting" only two quasiparticle states 
h = {BCS\ a_^a^{^ + ^ul,a+atj{um + Vmul,a+at^)-^ \BCS) 

Making an expansion to first order in a^ali^ of the inverse operator of the above formula (other terms containing 
more than two quasiparticle excitations are of course neglected) 

Replacing this quantity in I2 

h = {BCS\ a_^a^{^ + ^ul,a+alju-\l - VmUma+atm) \BCS) 
Again taking again into account only two quasiparticle excitations, 
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h = {BCS\ a_.^a^{ur.^ - Vm^)a+at„, \BCS) 



wc obtains: 

T — II t?iJni _ dUm 

knowing that the normalization condition is: 
we find by differentiation 

combining these two relations, we obtain in l2'. 

J _ 1 

-'2 ^ dpi 

then, the second term reads: 

h = (m, -ml (4)^^^^ \BCS) = -^f^ 

which can be cast as follows: 

\BCS)^--^6m (B2) 

The two matrix elements Ii and I2 corresponding to the two cases k 7^ I and k = I are now known. Reassembling the 
two parts Ii and I2 in the only one formula, we get: 

{M\ 4- \0) =h+l2 = (fc, -/| 4- \BCS) - (ukvi + uivk) {k\ ^ 10 (1 - 5kj) ~ ^^S,i 

Replacing this quantity in the cranking formula, noting that the crossed terms (/1/2 and I2I1) cancel in the product 
and knowing that Em — Eq = Ek + Ei (see subsec HTB]) we find: 

The first term of the r.h.s of the formula is the so-called non-diagonal term, whereas the second is the diagonal one. 
Appendix C: Another version of the cranking formula 

Let be H some nuclear hamiltonian and \0) , \M) its ground and excited state: noting that: 



(M| \0)^{Em-Eo){M\(^) ^ \0) 

\ / wave June \ / wave junc 

since the commutator gives: 

the cranking formula becomes: 

i/^i, ., /:^„| - 2fi J^^ (EM-Eof + ^f^^ [Em-Eo) 

where (i9/9/3i)u,ai,e func and (9/i9/3»)o^^„p p,,^^ have aheady defined in the appendix [B 

In the independent-particle approximation, we have in the second quantization representation: 

dH/dfii — ^ {v\ dHsp/dPi a^a^. After similar calculations than the ones performed in the appendix iBl we find 

h = (fc, -;| (9/9ft)»a.e func \BC S) = - [u^vi + (fc| dH^p/Op, \l) k^l 

therefore, the first part of the r.h.s of the above formula reads (the second remains unchanged): 

J^^ (EM-Eof - ^'^ (E,+E,r -WT 'WT ~ ^'''^> 

Then we obtain for the cranking formula: 

Appendix D: Final version of the cranking formula 

The expression 

E l 1 duk duk 
k 2Ek ^ dPi d0j 
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meet in the above formula can be further clarified. Recalling that: 

u, = (1/V2) (l+ek/V4+^y^' and V, = (l/V2) (l - Sk/V^+^Y^' 

where:e;j = — A is the single-particle energy with respect to the Fermi level and assuming for the moment that the 
deformation dependence appears through Cfc, A, and A, we have to evaluate ^ — f|^^ + ^0 + ^^ff • Thus, 
a simple differentiation of Uk with respect to leads to: 

> -1/2 

2V2 ^ 



} + [it (4 + ^r'^" - (4 + A^)"''' (-.it + At); 

multiplying by f^^and simplifying we get: 

aft - 2(4+A2) \^d0i ^kg^. J 

using Sk = Ck — A, we obtain: 

_ 1 / A^^f" - A - r^, _ \ 
^fc aft - 2(e2,+A2) X^Wi ^Wi ^^'^ '^''aft/ 

Moreover, noting that: 

it ~ 'wr l^-* 

we find: 

1 /2 

the quasiparticle energy is Ek = {e\ + A^) so that: 

-.-^lir = 2^- {a {k\ ^ |fc> - Af^ - (e, - A) |A} 
putting: 

- l^^l aft l'^) + Qp.+ A 
the product of the similar terms gives finally: 

r 1 1 duk 1 duk _ 1 A A 

^2 - Z^fe 2Efc 9ft «^ 9ft " l^k 2Ek ^ 2BJ^2SJ 

The cranking formula of the mass parameters becomes therefore: 

A, {A, M = ^ ^7"' +7/ {l\ ^ Ifc) (fcl ^ 10 (1 - 6,,) + 2fe J: §,RlR^ (Dl) 

[Ek + Ei) opi apj k k 
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